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Abstract. Motivated by the possible existence of other universes, with possible 
variations in the laws of physics, this paper explores the parameter space of 
fundamental constants that allows for the existence of stars. To make this problem 
tractable, we develop a semi-analytical stellar structure model that allows for physical 
understanding of these stars with unconventional parameters, as well as a means to 
survey the relevant parameter space. In this work, the most important quantities that 
determine stellar properties — and are allowed to vary — are the gravitational constant 
G, the fine structure constant a, and a composite parameter C that determines nuclear 
reaction rates. Working within this model, we delineate the portion of parameter 
space that allows for the existence of stars. Our main finding is that a sizable fraction 
of the parameter space (roughly one fourth) provides the values necessary for stellar 
objects to operate through sustained nuclear fusion. As a result, the set of parameters 
necessary to support stars are not particularly rare. In addition, we briefly consider 
the possibility that unconventional stars (e.g., black holes, dark matter stars) play the 
role filled by stars in our universe and constrain the allowed parameter space. 
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1. Introduction 

The current picture of inflationary cosmology allows for, and even predicts, the existence 
of an infinite number of space-time regions sometimes called pocket universes [H El 13] . I n 
many scenarios, these separate universes could potentially have different versions of the 
laws of physics, e.g., different values for the fundamental constants of nature. Motivated 
by this possibility, this paper considers the question of whether or not these hypothetical 
universes can support stars, i.e., long-lived hydrostatically supported stellar bodies that 
generate energy through (generalized) nuclear processes. Toward this end, this paper 
develops a simplified stellar model that allows for an exploration of stellar structure 
with different values of the fundamental parameters that determine stellar properties. 
We then use this model to delineate the parameter space that allows for the existence 
of stars. 

A great deal of previous work has considered the possibility of different values of the 
fundamental constants in alternate universes, or, in a related context, why the values of 
the constants have their observed values in our universe (e.g., [HE]). More recent papers 
have identified a large number of possible constants that could, in principle, vary from 
universe to universe. Different authors generally consider differing numbers of constants, 
however, with representative cases including 31 parameters [6] and 20 parameters [TJ. 
These papers generally adopt a global approach (see also [El El [TO]), in that they consider 
a wide variety of astronomical phenomena in these universes, including galaxy formation, 
star formation, stellar structure, and biology. This paper adopts a different approach 
by focusing on the particular issue of stars and stellar structure in alternate universes; 
this strategy allows for the question of the existence of stars to be considered in greater 
depth. 

Unlike many previous efforts, this paper constrains only the particular constants of 
nature that determine the characteristics of stars. Furthermore, as shown below, stellar 
structure depends on relatively few constants, some of them composite, rather than on 
large numbers of more fundamental parameters. More specifically, the most important 
quantities that directly determine stellar structure are the gravitational constant G, the 
fine structure constant a, and a composite parameter C that determines nuclear reaction 
rates. This latter parameter thus depends in a complicated manner on the strong and 
weak nuclear forces, as well as the particle masses. We thus perform our analysis in 
terms of this (a,G,C) parameter space. 

The goal of this work is thus relatively modest. Given the limited parameter space 
outlined above, this paper seeks to delineate the portions of it that allow for the existence 
of stars. In this context, stars are defined to be self-gravitating objects that are stable, 
long-lived, and actively generate energy through nuclear processes. Within the scope 
of this paper, however, we construct a more detailed model of stellar structure than 
those used in previous studies of alternate universes. On the other hand, we want 
to retain a (mostly) analytic model. Toward this end, we take the physical structure 
of the stars to be polytropes. This approach allows for stellar models of reasonable 
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accuracy; although it requires the numerical solution of the Lane-Emden equation, the 
numerically determined quantities can be written in terms of dimensionless parameters 
of order unity, so that one can obtain analytic expressions that show how the stellar 
properties depend on the input parameters of the problem. Given this stellar structure 
model, and the reduced (a, G, C) parameter space outlined above, finding the region of 
parameter space that allows for the existence of stars becomes a well-defined problem. 

As is well known, and as we re-derive below, both the minimum stellar mass and 
the maximum stellar mass have the same dependence on fundamental constants that 
carry dimensions [11] . More specifically, both the minimum and maximum mass can be 
written in terms of the fundamental stellar mass scale Mq defined according to 



where mp is the mass of the proton. As expected, the mass scale can be written as 
a dimensionless quantity {olq^ 2 ) times the proton mass; the appropriate value of the 
exponent (-3/2) in this relation is derived below. The mass scale Mq determines the 
allowed range of masses in any universe. 

In conventional star formation, our Galaxy (and others) produces stars with masses 
in the approximate range 0.08 < M*/M Q < 100, which corresponds to the range 
0.04 < M*/M < 50. One of the key questions of star formation theory is to understand, 
in detail, how and why galaxies produce a particular spectrum of stellar masses (the 
stellar initial mass function, or IMF) over this range [12]. Given the relative rarity of high 
mass stars, the vast majority of the stellar population lies within a factor of ~ 10 of the 
fundamental mass scale Mq. For completeness we note that the star formation process 
does not involve thermonuclear fusion, so that the mass scale of the hydrogen burning 
limit (at 0.08 M Q ) does not enter into the process. As a result, many objects with 
somewhat smaller masses - brown dwarfs - are also produced. One of the objectives 
of this paper is to understand how the range of possible stellar masses changes with 
differing values of the fundamental constants of nature. 

This paper is organized as follows. We construct a polytropic model for stellar 
structure in §2, and identify the relevant input parameters that determine stellar 
characteristics. Working within this stellar model, we constrain the values of the stellar 
input parameters in §3; in particular, we delineate the portion of parameter space that 
allow for the existence of stars. Even in universes that do not support conventional 
stars, those generating energy via nuclear fusion, it remain possible for unconventional 
stars to play the same role. These objects are briefly considered in §4 and include black 
holes, dark matter stars, and degenerate baryonic stars that generate energy via dark 
matter capture and annihilation. Finally, we conclude in §5 with a summary of our 
results and a discussion of its limitations, including an outline for possible future work. 




(1) 



where olq is the gravitational fine structure constant 
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2. Stellar Structure Models 

In general, the construction of stellar structure models requires the specification 
and solution of four coupled differential equations, i.e., force balance (hydrostatic 
equilibrium), conservation of mass, heat transport, and energy generation. This set 
of equations is augmented by an equation of state, the form of the stellar opacity, 
and the nuclear reaction rates. In this section we construct a polytropic model of stellar 
structure. The goal is to make the model detailed enough to capture the essential physics 
and simple enough to allow (mostly) analytic results, which in turn show how different 
values of the fundamental constants affect the results. Throughout this treatment, we 
will begin with standard results from stellar structure theory [Til D21 EE] and generalize 
to allow for different stellar input parameters. 



2.1. Hydrostatic Equilibrium Structures 

In this case, we will use a polytropic equation of state and thereby replace the force 
balance and mass conservation equations with the Lane-Emden equation. The equation 
of state thus takes the form 

P = Kp r where T = 1 + - , (3) 

n 

where the second equation defines the polytropic index n. Note that low mass stars 
and degenerate stars have polytropic index n = 3/2, whereas high mass stars, with 
substantial radiation pressure in their interiors, have index n — > 3. As a result, the 
index is slowly varying over the range of possible stellar masses. Following standard 
methods pH QH Q31 Q3] , we define 

r KT 

S = p=p c f, and R = (r _ 1)4;rG) , c2 _ r . (4) 

so that the dimensionless equation for the hydrostatic structure of the star becomes 
d_ (pdf 

Here, the parameter p c is the central density (in physical units) so that / n (£) is the 
dimensionless density distribution. For a given polytropic index n (or a given T), 
equation §5§ thus specifies the density profile up to the constants p c and R. Note that 
once the density is determined, the pressure is specified via the equation of state 
Further, in the stellar regime, the star obeys the ideal gas law so that the temperature 
is given by T = P/(JZp), with 1Z = k/(m); the function /(£) thus represents the 
dimensionless temperature profile of the star. Integration of equation (jSJ) outwards, 
subject to the boundary conditions / = 1 and df/d£ = at £ = 0, then determines 
the position of the outer boundary of the star, i.e., the value where /(£*) = 0. As a 
result, the stellar radius is given by 

R* = R^ . (6) 



+er = o. (5) 
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Figure 1. Density, pressure, and temperature distributions for n = 3/2 polytrope. 
The solid curve shows the density profile p{£)/p c , the dashed curve shows the pressure 
profile P(£)/P c , and the dotted curve shows the temperature profile /(£) = T(£)/T c . 
For a polytrope, the variables are related through the expressions P oc p 1+1 / n and 
pcx/". 

The physical structure of the star is thus specified up to the constants p c and R. 
These parameters are not independent for a given stellar mass; instead, they are related 
via the constraint 

M* = 4ttR 3 Pc £ 2 / n (£R = 47ri? 3 P c^o , (7) 
Jo 

where the final equality defines the dimensionless quantity fiQ, which is of order unity 
and depends only on the polytropic index n. 
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2.2. Nuclear Reactions 

The next step is to estimate how the nuclear ignition temperature depends on more 
fundamental parameters of physics. Thermonuclear fusion generally depends on three 
physical variables: the temperature T, the Gamow energy Eq, and the nuclear fusion 
factor S(E). The Gamow energy is given by 

E G = (naZ^) 2 2mim2 c 2 = ( 7ra Z 1 Z 2 ) 2 2m R c 2 , (8) 
mi + ni2 

where rrij are the masses of the nuclei, Zj are their charge (in units of e), and where the 
second equality defines the reduced mass. For the case of two protons, E G = 493 keV. 
The parameter a is the usual (electromagnetic) fine structure constant 

a = Yc~m> (9) 

where the numerical value applies to our universe. Thus, the Gamow energy, which 
sets the degree of Coulomb barrier penetration, is determined by the strength of the 
electromagnetic force (through a). The strength of the strong and weak nuclear forces 
enter into the problem by setting the nuclear fusion factor S(E), which in turn sets the 
interaction cross section according to 

■E G \^ 



(F\ S W 
a{E) = — ^- exp 



■e 



(10) 



where E is the energy of the interacting nuclei. The temperature at the center of the 
star determines the distribution of E. Under most circumstances in ordinary stars, the 
cross section has the approximate dependence a oc 1/E so that the nuclear fusion 
factor S(E) is a slowly varying function of energy. This dependence arises when 
the cross section is proportional to the square of the de Broglie wavelength, so that 
a ~ A 2 ~ {h/p) 2 ~ h 2 /(2mE); this relation holds when the nuclei are in the realm of 
non-relativistic quantum mechanics. 

The nuclei generally have a thermal distribution of energy so that 

/ 8 \ 1/2 / 1 \ 3/2 r°° 

<™> = bj fe) /„ °(e)°M-Elk T] EAE. (11) 

As a result, the effectiveness of nuclear reactions is controlled by an exponential factor 
exp[— $], where the function $ has contributions from the cross section and the thermal 
distribution, i.e., 

* = + • < 12) 

The integral in equation (TTTj) is dominated by energies near the minimum of $, where 
E = Eq = EQ 3 (kT/2) 2 ^ 3 , and where the function takes the value 

/ F \ 1 / 3 
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If we approximate the integral using Laplace's method [16], the reaction rate R\ 2 for 
two nuclear species with number densities n\ and n 2 can be written in the form 

R12 = nm 2 8 S(E )e 2 exp[-3Q] , (14) 

Vo7ro;ZiZ 2 m j RC 

where we have defined 

MIT- 

2.3. Stellar Luminosity and Energy Transport 

The luminosity of the star is determined through the equation 

Jl = Arrr 2 e(r) , (16) 
dr 

e is the luminosity density, i.e., the power generated per unit volume. This quantity can 
be written in terms of the nuclear reaction rates via 

e{r) =Cp 2 e 2 exp[-36], (17) 

where is defined above, and where 

c = (AE)R 12 exp[3@] = 8(AE)S(E ) 

p 2 6 2 ^/3irami , m 2 Z 1 Z 2 mRC ' 

where (AE) is the mean energy generated per nuclear reaction. In our universe 
C ~ 2 x 10 4 cm 5 s -3 g^ 1 for proton-proton fusion under typical stellar conditions. 
The total stellar luminosity is given by the integral 

L* = C4nR 3 p c 2 [^ f n £ 2 e 2 exp[-3QW = C^R 3 p 2 I{Q c ) , (19) 
Jo 

where the second equality defines I(O c ), and where 9 C = 0(£ = 0) = (Eq /AkTc) 1 ' 3 . Note 
that for a given polytrope, the integral is specified up to the constant C : T = T c f(£), 

e = e c /-V3 ( e). 

At this point, the definition of equation (j4]), the mass integral constraint ([7j), and 
the luminosity integral ( JT9l) provide us with three equations for four unknowns: the 
radial scale R, the central density p c , the total luminosity L*, and the coefficient K in 
the equation of state. Notice that if the star is degenerate, then the coefficient K is 
specified by quantum mechanics, V = 5/3, and one could solve the first two of these 
equations for R and p c , thereby determining the physical structure of the star. Note 
that the quantum mechanical value of K represents the minimum possible value. If the 
star is not degenerate, but rather obeys the ideal gas law, then the central temperature 
is related to the central density through 1ZT C = Kp c l l n , so that T c does not represent 
a new unknown, and the stellar luminosity L* is the only new variable introduced by 
luminosity equation (fl9j) . 

For ordinary stars, one needs to use the fourth equation of stellar structure to finish 
the calculation. In the case of radiative stars, the energy transport equation takes the 
form 

^ndT 2>pn L(r) , , 

dr Aac 4nr z 
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L* L -^d^^^—R, (21) 



where k is the opacity. In the spirit of this paper, we want to obtain a simplified set 
of stellar structure models to consider the effects of varying constants. As a result, we 
make the following approximation. The opacity k generally follows Kramer's law so that 
pT~ 7 / 2 . For the case of polytropic equations of state, we find that up ~ p . 
For the particular case n = 7/4, the product up is strictly constant. For other values 
of the polytropic index, the quantity np is slowly varying. As a result, we assume up = 
Kop c = constant for purposes of solving the energy transport equation (120]) . This ansatz 
implies that 

where we have defined = L(£)/L*. The full expression for is given by the 
integral in equation (fl9|) . For purposes of solving equation (12~TI) . however, we make a 
further simplification: We assume that the integrand of equation ( fl9l) is sharply peaked 
toward the center of the star, and that the nuclear reaction rates depend on a power- 
law function of temperature. Consistency then demands that the power-law index is C . 
Further, the temperature can be modeled as an exponentially decaying function near 
the center of the star so that T ~ exp [—/?£]. The expression for £(£) then becomes 

= o r^^e^dx where x end = /39 c £ . (22) 
2 jo 

Using this expression for £(£) in the integral of equation ( 12~TT) . we can write the luminosity 
in the form 

A 47TC R , . 

U = aT c \ — . 23 

2.4- Stellar Structure Solutions 

With the solution (1231) to the energy transport equation, we now have four equations 
and four unknowns. After some algebra, we obtain the following equation for the central 
temperature 

or, alternately, 

2 12 7T 5 

45 (3 Ko CE^ti 3 

The right hand side of the equation is thus a dimensionless quantity. Further, the 
quantities fiQ and (3 are dimensionless measures of the mass and luminosity integrals 
over the star, respectively; they are expected to be of order unity and to be roughly 
constant from star to star (and from universe to universe). The remaining constants are 
fundamental. Note that for typical values of the parameters in our universe, the right 
hand side of this equation is approximately 10 -9 . 




/(e c )e c - 8 = AC - _,, 3 2 1 — 1 (z 1 ^) • (25) 
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With the central temperature T c , or equivalently, 9 C , determined through equation 
(j25p . we can find expressions for the remaining stellar parameters. The radius is given 
by 

GM^m) 

kT c (n + l)/i 
and the luminosity is given by 

_ 167T 4 1 (MA 3 fG( m )V 

15 tfc 2 (3K Q c \no) \n + l) 1 J 

The photospheric temperature is then determined from the usual outer boundary 
condition so that 

M^J • (28) 

For this simple polytropic stellar model, Figures [2] and [3] show the H-R diagram 
and the corresponding luminosity versus mass relation for stars on the zero age main 
sequence (ZAMS). The three curves show different choices for the polytropic indices: 
The dashed curves show results for n = 3/2, the value appropriate for low-mass stars. 
The dotted curves show the results for n = 3, the value for high-mass stars. The solid 
line (marked by symbols) show the results for n varying smoothly between n = 3/2 
in the limit M* — > and n = 3 in the limit M* — > oo. We take this latter case as 
our standard model (although the effects of changing the polytropic index n are small 
compared to the effects of changing the fundamental constants - see §3). 

One can compare these models with the results of more sophisticated stellar 
structure models ( [T3l [13] ) or with observations of stars on the ZAMS. In both of these 
comparison, this polytropic model provides a good prediction for the stellar temperature 
as a function of stellar mass. However, the luminosities of the highest mass stars are 
somewhat low, mostly because the stellar radii from the models are correspondingly 
low; this discrepancy, in turn, results from our simplified treatment of nuclear reactions. 
Nonetheless, this polytropic model works rather well, and produces the correct stellar 
characteristics (L*, i?*, T#), within a factor of ~ 2, as a function of mass M*, over a 
range in mass of ~ 1000 and a range in luminosity of ~ 10 10 . This degree of accuracy 
is sufficient for the purposes of this paper, and is quite good given the simplifying 
assumptions used in order to obtain analytic results. More sophisticated stellar models 
would include varying values of C to incorporate more complex nuclear reaction chains, 
detailed energy transport including convection, a more refined treatment of opacity, 
and a fully self-consistent determination of the density and pressure profiles (i.e., the 
departures from our polytropic models). In particular, we can achieve even better 
agreement between this stellar structure model and observed stellar properties if we 
allow the nuclear reaction parameter C to increase with stellar mass (as it does in high 
mass stars due to the CNO cycle). In the spirit of this work, however, we use a single 
value of C, which corresponds to the case in which a single nuclear species is available 
for fusion (this scenario thus represents the simplest universes). 



Stars In Other Universes: Stellar structure with different fundamental constants 10 




Figure 2. H-R Diagram showing the main sequence for polytropic stellar model using 
standard values of the parameters, i.e., those in our universe. The three cases shown 
here correspond to the main sequence for an n = 3/2 polytrope (lower dashed curve), 
an n = 3 polytrope (upper dotted curve), and a model that smoothly varies from n = 
3/2 at low masses to n — 3 at high masses (solid curve marked by symbols). 



3. Constraints on the Existence of Stars 

Using the stellar structure model developed in the previous section, we now explore 
the range of possible stellar masses in universes with varying value of the stellar 
parameters. First, we find the minimum stellar mass required for a star to overcome 
quantum mechanical degeneracy pressure (§3.1) and then find the maximum stellar 
mass as limited by radiation pressure (§3.2). These two limits are then combined to 
find the allowed range of stellar masses, which can vanish when the required nuclear 
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Figure 3. Stellar luminosity as a function of stellar mass for standard values of the 
parameters. The three curves shown here correspond to the — relation for an 
n = 3/2 polytrope (dashed curve), an n = 3 polytrope (dotted curve), and a model 
that smoothly varies from n = 3/2 at low masses to n = 3 at high masses (solid curve 
marked by symbols). All quantities are given in solar units. 



burning temperatures becomes too high (§3.3). Another constraint on stellar parameters 
arises from the requirement that stable nuclear burning configurations exist (§3.4). We 
delineate (in §3.5) the range of parameters for which these two considerations provide 
the limiting constraints on stellar masses and then find the region of parameter space 
that allows the existence of stars. Finally, we consider the constraints implied by the 
Eddington luminosity (§3.6) and show that they are comparable to those considered in 
the previous subsections. 
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3.1. Minimum Stellar Mass 

The minimum mass of a star is determined by the onset of degeneracy pressure. 
Specifically, for stars with sufficiently small masses, degeneracy pressure enforces a 
maximum temperature which is below that required for nuclear fusion. The central 
pressure at the center of a star is given approximately by the expression 

Pc-(^) 1/3 GM*/V /3 , (29) 

where the subscript denotes that the quantities are to be evaluated at the center of the 
star. This result follows directly from the requirement of hydrostatic equilibrium (e.g., 
[15]). 

At the low mass end of the range of possible stellar masses, the pressure is 
determined by contributions from the ideal gas law and from non-relativistic electron 
degeneracy pressure. As a result, the central pressure of the star must also satisfy the 
relation 

-^-)kT c + K d J^-) 5/3 , (30) 

where m ion is the mean mass of the ions (so that p c /m ion determines the number density 
of ions) and where the constant K dp that determines degeneracy pressure is given by 

h 2 (n 2 ^2/3 



A' dp = — (3^j , (31) 

where m e is the electron mass. Notice that we have also assumed that the star has 
neutral charge so that the number density of electrons is equal to that of the ions, and 
that m e <C m ion . 

Combining the two expressions for the central pressure and solving for the central 
temperature, we obtain 

1 /3 

kT c = (^j GMl' z m ionPc l ' z - K dp ( Pc /m ion ) 2 / 3 . (32) 

The above expression is a simple quadratic function of the variable Pc 1 ^ 3 and has a 
maximum for a particular value of the central density i.e., 

kT max = ( — ) - M * miOD . (33) 
V36/ 4K dp v ; 

If we set this value of the central temperature equal to the minimum required ignition 
temperature for a star, T nuc , we obtain the minimum stellar mass 



,3/4 



36\ 1/2 (4A' dp 4T, 



"ion 



After rewriting the equation of state parameter K dp in terms of fundamental constants, 
this expression for the minimum stellar mass becomes 
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As expected, the minimum stellar mass is given by a dimensionless expression times the 
fundamental stellar mass scale defined in equation ([I]). Notice also that the gravitational 
constant G enters into this mass expression with an exponent of -3/2, as anticipated by 
equation ([!]). 



3.2. Maximum Stellar Mass 

A similar calculation gives the maximum possible stellar mass. In this case the central 
pressure also has two contributions, this time from the ideal gas law and from radiation 
pressure Pr, where 

P R =\aT c \ (36) 

where a = ir 2 k 4 /15(hc) 3 is the radiation constant. Following standard convention [11], 
we define the parameter f g to be the fraction of the central pressure provided by the ideal 
gas law. As a result, the radiation pressure contribution is given by Pr = (1 — f g )P c - 
The central temperature can be eliminated in favor of f g to obtain the expression 



!. ) VW. 

where (m) is the mean mass per particle of a massive star. By demanding that the star 
is in hydrostatic equilibrium, we obtain the following expression for the maximum mass 
of a star: 

^-(f)"(§<^)^(£)\ 

which can also be written in terms of the fundamental mass scale M , i.e., 



where this expression must be evaluated at the maximum value of f g for which the 
star can remain stable. Although the requirement of stability does not provide a 
perfectly well-defined threshold for f g , the value f g = 1/2 is generally used [11] and 
predicts maximum stellar masses in reasonable agreement with observed stellar masses 
(for present-day stars in our universe). For this choice, the above expression becomes 
M, max « 20(mp/(m)) 2 M . Since massive stars are highly ionized, (m) pa O.Gmp under 
standard conditions, and hence M* max ps 56M ~ 100M© for our universe. As shown 
below, this constraint is nearly the same as that derived on the basis of the Eddington 
luminosity (§3.6). 



3.3. Constraints on the Range of Stellar Masses: 
The Maximum Nuclear Ignition Temperature 

As derived above, the minimum stellar mass can be written as a dimensionless coefficient 
times the fundamental stellar mass scale from equation ([1]). Further, the dimensionless 
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coefficient depends on the ratio of the nuclear ignition temperature to the electron mass 
energy, i.e., kT nnc /m e c 2 . The maximum stellar mass, also defined above, can be written 
as a second dimensionless coefficient times the mass scale Mq. This second coefficient 
depends on the maximum radiation pressure fraction / g and (somewhat less sensitively) 
on the mean particle mass (m) of a high mass star. For completeness, we note that the 
Chandrasekhar mass M c h [15J can be written as yet another dimensionless coefficient 
times this fundamental mass scale, i.e., 

M ch « \ (2vr) 3 / 2 (I) 2 Mo, (40) 

where Z/A specifies the number of electrons per nucleon in the star. 

These results thus show that if the constants of the universe were different, or if 
they are different in other universes (or different in other parts of our universe), then 
the possible range of stellar masses would change accordingly. We see immediately that 
if the nuclear ignition temperature is too large, then the range of stellar masses could 
vanish. If all other constants are held fixed, then the requirement that the minimum 
stellar mass becomes as large as the maximum stellar mass is given by 

4/3 



( kT nuc ^ 




/360\ 






UttV 




\ 8/3 

rru, ■ 



- > , \ - \ i « 1-4 > ( 41 ) 

V8f g 2 ) \ (m) ) \{m)) 1 ; 

where we have used / g = 1/2 to obtain the final equality. For high mass stars in our 
universe, {m)/m- lon = 0.6, and the right hand side of the equation is about 5.6. For the 
simplistic case where (m) = m = m ion , the right hand side is 1.4. In any case, this value 
is of order unity and is not expected to vary substantially from universe to universe. As 
a result, the condition for the nuclear burning temperature to be so high that no viable 
range of stellar masses exists takes the form kT nuc / (m e c 2 ) ^ 2. For standard values of the 
other parameters, the nuclear ignition temperature (for Hydrogen fusion) would have 
to exceed T nuc ~ 10 10 K. For comparison, the usual Hydrogen burning temperature 
is about 10 7 K and the Helium burning temperature is about 2 x 10 s K. We stress 
that the Hydrogen burning temperature in our universe is much smaller than the value 
required for no range of stellar masses to exist — in this sense, our universe is not fine- 
tuned to have special values of the constants to allow the existence of stars. The large 
value of nuclear ignition temperature required to suppress the existence of stars roughly 
corresponds to the temperature required for Silicon burning in massive stars (again, for 
the standard values of the other parameters). Finally we note that the nuclear burning 
temperature T nuc depends on the fundamental constants in a complicated manner; this 
issue is addressed below. 

Equation (141j) emphasizes several important issues. First we note that the existence 
of a viable range of stellar masses — according to this constraint — does not depend 
on the gravitational constant G. The value of G determines the scale for the stellar 
mass range, and the scale is proportional to G~ 3 ^ 2 ~ olq^" 1 \ but the coefficients that 
define both the minimum stellar mass and the maximum stellar mass are independent of 
G. The possible existence of stars in a given universe depends on having a low enough 
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nuclear ignition temperature, which requires the strong nuclear force to be "strong 
enough" and/or the electromagnetic force to be "weak enough". These requirements 
are taken up in §3.5. Notice also that we have assumed m e C mp, so that electrons 
provide the degeneracy pressure, but the ions provide the mass. 



3.4- Constraints on Stable Stellar Configurations 

In this section we combine the results derived above to determine the minimum 
temperature required for a star to operate through the burning of nuclear fuel (for 
given values of the constants). For a given minimum nuclear burning temperature 
^nuo equation (|35|) defines the minimum mass necessary for fusion. Alternatively, the 
equation gives the maximum temperature that can be attained with a star of a given 
mass in the face of degeneracy pressure. On the other hand, equation (125]) specifies the 
central temperature T c necessary for a star to operate as a function of stellar mass. We 
also note that the temperature T c is an increasing function of stellar mass. By using 
the minimum mass from equation (1351) to specify the mass in equation (j25l) . we can 
eliminate the mass dependence and solve for the minimum value of the nuclear ignition 
temperature T nuc . The resulting temperature is given in terms of G c , which is given by 
the solution to the following equation: 

— » - m © fe) fey s ■ 

Note that the parameters on the right hand side of the equation have been grouped to 
include numbers, constants that set units, dimensionless parameters of the polytropic 
solution, the relevant particle masses, and the stellar parameters that depend on the 
fundamental forces. Within the treatment of this paper, these latter quantities could 
vary from universe to universe. Notice also that we have specialized to the case in which 
(m) = m ion = m. 

The left hand side of equation (I42p is determined for a given polytropic index. 
Here we use the value n = 3/2 corresponding to both low- mass conventional stars and 
degenerate stars. The resulting profile for C /(0 C ) is shown in Figure HI The right hand 
side of equation (j4"2"l) depends on the fundamental constants and is thus specified for a 
given universe. In order for nuclear burning to take place, equation (|4"2l) must have a 
solution — the left hand side has a maximum value, which places an upper bound on 
the parameters of the right hand side. Through numerical evaluation, we find that this 
maximum value is ~0.0478 and occurs at @ c ~ 0.869. The maximum possible nuclear 
burning temperature thus takes the form 

(kT) max « 0.38£ G , (43) 

where Eq is the Gamow energy appropriate for the given universe. The corresponding 
constraint on the stellar parameters required for nuclear burning can then be written in 
the form 



tfG 5 n /?/i 4 



c 2 mmlK C 2 23 7r 7 3 4 



e c /(e c ) 



2.6 x 10" 5 , (44) 
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where we have combined all dimensionless quantities on the right hand side. For typical 
stellar parameters in our universe, the left hand side of the above equation has the 
value ~ 2.4 x 10~ 9 , smaller than the maximum by a factor of ~1 1,000. As a result, 
the combination of constants derived here can take on a wide range of values and still 
allow for the existence of nuclear burning stars. In this sense, the presence of stars in 
our universe does not require fine-tuning the constants. 

Notice that for combinations of the constants that allow for nuclear burning, 
equation f j4"2l) has two solutions. The relevant physical solution is the one with larger 
G c , which corresponds to lower temperature. The second, high temperature solution 
would lead to an unstable stellar configuration. As a consistency check, note that for 
the values of the constants in our universe, the solution to equation fj4"2l implies that 
G c ~ 5.38, which corresponds to a temperature of about 9 x 10 6 K. This value is thus 
approximately correct: Detailed stellar models show that the central temperature of the 
Sun is about 15 x 10 6 K, and the lowest possible hydrogen burning temperature is a few 
million degrees [TTJ, [131 E! • 

3.5. Combining the Constraints 

Thus far, we have derived two constraints on the range of stellar structure parameters 
that allow for the existence of stars. The requirement of stable nuclear burning 
configuration places an upper limit on the nuclear burning temperature, which takes 
the approximate form kT ;$ 0.38Eq. In addition, the requirement that the minimum 
stellar mass (due to degeneracy pressure) not exceed the maximum stellar mass (due 
to radiation pressure) places a second upper limit on the nuclear burning temperature, 
kT £ 2m e c 2 . As a result, the reason for a universe failing to produce stars depends on 
the size of the dimensionless parameter 

Q F = 2a 2 — , (45) 

m e 

where m is the mass of the nuclei that would experience reactions. Note that Qp is 
proportional to the ratio of the Gamow energy to the rest mass energy of the electron 
and has the value Qp ~ 0.2 in our universe. For Qp > 1, stars can fail to exist due to 
the range of allowed stellar masses shrinking to zero, whereas for Qp < 1 stars can fail 
to exist due to the absence of stable nuclear burning configurations. 

We can combine the constraints to delineate the portion of parameter space that 
allows for the existence of stars. For the sake of definiteness, we fix the values of the 
particle masses and specialize to the simplest case where the nuclear burning species has 
a single mass m. We also assume that the stellar opacity scales according to k oc a 2 , 
as expected since k ~ apjm and o~p oc a 2 . With these restrictions, the remaining 
stellar parameters that can be varied are the fine structure constant a, the gravitational 
constant G, and the nuclear burning parameter C. Note that a depends on the strength 
of the electromagnetic force, G depends on the strength of gravity, and C depends on a 
combination of the weak and strong nuclear forces, which jointly determine the nuclear 
reaction properties for a given universe. Notice also that since we are fixing particle 
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Figure 4. Profile of e c 7(6 c ) as a function of 6 C = (£ , G /4fcT c ) 1 / 3 . The integral J(6 C ) 
determines the stellar luminosity in dimensionless units and O c defines the central 
stellar temperature. This profile has a well-defined maximum near C ss 0.869, where 
the peak of the profile defines a limit on the values of the fundamental constants 
required for nuclear burning, and where the location of the peak defines a maximum 
nuclear burning temperature (see text). 

masses, the gravitational constant G is proportional to the gravitational fine structure 
constant otc (eq. [2]). 

Figure [5] shows the resulting allowed region of parameter space for the existence of 
stars. Here we are working in the (a, G) plane, where we scale the parameters by their 
values in our universe, and the results are presented on a logarithmic scale. For a given 
nuclear burning constant C, Figure [5] shows the portion of the plane that allows for stars 
to successfully achieve sustained nuclear reactions. Curves are given for three values of 
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Figure 5. Allowed region of parameter space for the existence of stars. Here the 
parameter space is the plane of the gravitational constant log 10 [G/Go] versus the fine 
structure constant log 10 [a:/ao], where both quantities are scaled relative to the values 
in our universe. The allowed region lies under the curves, which are plotted here for 
three different values of the nuclear burning constants C: the standard value for p-p 
burning in our universe (solid curve), 100 times the standard value (dashed curve), 
and 0.01 times the standard value (dotted curve). The open triangular symbol marks 
the location of our universe in this parameter space. 
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C: the value for p-p burning in our universe (solid curve), 100 times larger than this 
value (dashed curve), and 100 times smaller (dotted curve). The region of the diagram 
that allows for the existence of stars is the area below the curves. 

Figure [5] provides an assessment of how "fine-tuned" the stellar parameters must 
be in order to support the existence of stars. First we note that our universe, with 
its location in this parameter space marked by the open triangle, does not lie near the 
boundary between universes with stars and those without. Specifically, the values of a, 
G, and/or C can change by more than two orders of magnitude in any direction (and 
by larger factors in some directions) and still allow for stars to function. This finding 
can be stated another way: Within the parameter space shown, which spans 10 orders 
of magnitude in both a and G, about one fourth of the space supports the existence of 
stars. 

Next we note that a relatively sharp boundary occurs in this parameter space for 
large values of the fine structure constant, where a ~ 200a , and this boundary is nearly 
independent of the nuclear burning constant C. Strictly speaking, this well-defined 
boundary is the result of the required value of G becoming an exponentially decreasing 
function of a/atQ, as shown in §3.7 below. For the given range of G and for values of 
a above this threshold, the Gamow energy is much larger than the rest mass energy 
of the electron, so that the maximum nuclear burning temperature becomes a fixed 
value (that given by eq. [41]). and hence the nuclear reaction rates are exponentially 
suppressed by the electromagnetic barrier (§2.2). On the other side of the graph, for 
values of a smaller than those in our universe, the range of allowed parameter space 
is limited due to the absence of stable nuclear burning configurations (§3.4). In this 
regime, for sufficiently large G, the nuclear burning temperature becomes so large that 
the barrier disappears (and hence stability is no longer possible). Since the nuclear 
burning temperature T nuc required to support stars against gravity increases as the 
gravitational constant G increases, and since T nuc is bounded from above, there is a 
maximum value of G that can support stars (for a given value of C). For the value of C 
appropriate for p-p burning in our universe, we thus find that G/Gq £ 2 x 10 5 . Finally, 
we note that "stellar" bodies outside the range of allowed parameter space can exist, 
in principle, and can even generate energy, but they would not resemble the stable, 
long-lived nuclear burning stars of our universe. 

3.6. The Eddington Luminosity 

For a star of given mass, the maximum rate at which it can generate energy is given by 
the Eddington luminosity. This luminosity defines a minimum lifetime for stars. The 
Eddington luminosity can be written in the form 

£*max = ^CGM^/ K em , (46) 

where K em is the opacity in the stellar photosphere. For the sake of definiteness, we take 
K em to be the opacity appropriate for pure electron scattering, which is applicable to 
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hot plasmas where the Eddington luminosity is relevant, i.e., 
1+Xi a T 



(47) 



"till f-. 7 

2 mp 

where o~t is the Thompson cross section and Xi is the mass fraction of Hydrogen. Since 
the maximum luminosity implies a minimum stellar lifetime, for a given efficiency e of 
converting mass into energy, we obtain the following constraint on stellar lifetimes 

/l + Xi\ a 2 hmp 

U > **min = e — — ^ . 48 

V 3 / (xq {m e c) z 

Since atomic time scales are given (approximately) by 

a z m e c z 

the ratio of stellar time scales to atomic time scales is given by the following expression: 

^*min 

/l + X±\ a 4 m P 



\ > / ( 5 °) 

t A \ 6 J a G m e 
where the expression has a numerical value of ~ 4 x 10 30 for the parameters in our 
universe. 

We can also use the Eddington luminosity to derive another upper limit on the 
allowed stellar mass. Within the context of our model, the stellar luminosity is given 
by equation (1271) . This luminosity must be less than the Eddington luminosity given by 
equation ( 14*61) . which implies a constraint of the form 

M„ < 4 ^ / [3fi 3 K m P O c \ 1/2 ^ (51) 



Mo ~ 7T V 0~T 

where we have specialized to the case of polytropic index n = 3 (appropriate for high 

mass stars with large admixtures of radiation pressure) and have taken (m) = mp. Note 

that since k ~ ctt/^p, and since (3 and /^o are given by the polytropic solution (and 

are of order unity), the right hand side of the above equation is approximately 5O\/0c, 

as expected. In other words, the requirement that the stellar luminosity must be less 

than the Eddington limit (eq. [46]) produces nearly the same bound on stellar masses 

as the requirement that the star not be dominated by radiation pressure (eq. [39]). 

Notice also that we expect k,q ~ ax/iTLp for other universes, so that the general 

constraint takes the approximate form M*/M £ 50 a/©!- In addition, as shown by 

Figure HJ the parameter C is confined to a narrow range — the function B C J(G C ), 

and hence the left hand side of equation (I42p . varies by 8 orders of magnitude for 
i < .far < q 



3. 7. Limiting Forms 

For much of the allowed parameter space where stars can operate; the value of C is 
large compared to its minimum value; specifically, this claim holds for the region of 
parameter space that is not near the upper left boundary in Figure [5j In this case, 
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we can derive an analytic asymptotic expression for the integral function J(0 C ), which 
takes the form 



Comparing this asymptotic expression to the numerically determined values, we find 
that equation (152]) provides an estimate that is within a factor of 2 of the correct result 
over the range 1 < C < 100, where /(O c ) varies by ~ 128 orders of magnitude. 

With this asymptotic expression in hand, we can find the relationship between the 
gravitational constant and the fine structure constant on the boundary of parameter 
space (as shown in Figure [5]). We find that 



At the edge of the allowed stellar parameter space, G is thus an exponentially decreasing 
function a, which results in the nearly vertical boundary shown in Figure O 

4. Unconventional Stars 

The results of the previous section show that stars can exist in a relatively large fraction 
of the parameter space. On the other hand, in order for stars to exist at all, the nuclear 
burning parameter C must be nonzero; otherwise, stars, as objects powered by nuclear 
reactions, cannot exist. In situations where C = 0, or where the values of the other 
parameters conspire to disallow stars (see Fig. [5]), other types of stellar objects could, 
in principle, fill the role played by stars in our universe. This section briefly explores 
this possibility with three examples: black holes that generate energy through Hawking 
evaporation (§4.1), degenerate dark matter stars that generate energy via annihilation 
(§4.2), and degenerate baryonic matter stars that generate energy by capturing dark 
matter particles which then annihilate (§4.3). We note that a host of other possibilities 
exist (e.g., astrophysical objects powered by proton decay), but a proper treatment of 
such cases is beyond the scope of this present work. 

4-1. Black Holes 

Black holes can exist in any universe with gravity and will generate energy (at some 
rate) through Hawking evaporation (e.g., [17] )■ Further, the stellar structure of these 
objects depends only on the gravitational constant G. In order to consider black holes 
playing the role of stars, however, we must invoke additional constraints. For purposes 
of illustration, this section finds the values of the fundamental constants for which black 
holes can serve as stellar bodies to support "life" . Specifically, in order for black holes 
to fill the role played by stars in our universe, two constraints must be satisfied: First, 
the black holes must live long enough to allow for life to develop. Second, the black 
holes must provide enough power to run a biosphere. The first constraint implies that 
black holes must be sufficiently massive, whereas the second constraint implies that the 




(52) 




(53) 
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black holes must be sufficiently small. The compromise between these two requirements 
provides an overall constraint that must be met in order for black holes to play the role 
of stars. 

The lifetime of a black hole with mass M b h is given by 

26507T ^n,^ o . . 

nn = -^~ T G 2 M hh 3 , (54) 

where g* is the total number of effective degrees of freedom in the radiation field 
produced through the Hawking effect. This lifetime should be compared with the typical 
atomic time scale ta given by equation (jUJ). We thus have a constraint of the form 

— = ^^« 2 G 2 M bh 3 m e > N hio , (55) 

where iV bio is the number of atomic time scales required for life to evolve. In our solar 
system, the number N^ ~ 10 34 , which is also the number of atomic time scales in 
the life of a solar type star. Although the minimum value of iV biC) remains uncertain, 
we expect it to be within a few orders of magnitude of this value. For the sake of 
definiteness, we take the (somewhat optimistic) value of A^ bio = 10 33 for this analysis. 

The second constraint is that the black hole must provide enough power to run 
a biosphere. In our solar system, the Earth intercepts about 100 quadrillion Watts 
of power from the Sun. We thus expect that the black hole must have a minimum 
luminosity and obey a constraint of the form 

L hh = ^f(GM hh y 2 >L min , (56) 
76807T 

where L min is the minimum luminosity of a stellar object required to support life. In 
general, this minimum value of luminosity will vary with the values of the fundamental 
constants. In the absence of a definitive theory, we adopt the following simple scaling 
law: The energy levels Ea of atoms vary according to Ea oc a 2 , and the atomic time 
scale varies as t& oc a~ 2 . In order for the luminosity to provide the same number of 
atomic reactions over the total lifetime of the system, the luminosity should scale with 
the fine structure constant as 

£min = L min0 (a/a ) A , (57) 

where L m i n0 is the minimum necessary luminosity in our universe. Although the value 
of this latter quantity is uncertain, we adopt L mm0 « 10 17 erg s _1 as a representative 
value. The scaling law of equation ( 1571) is also not definitive, but rather illustrative. 

Combining the two constraints allows for the elimination of the mass, and thereby 
provides an overall constraint of the from 



N hio h 1/2 (G/a 



,2 „7 



< ; 57? • (58) 

m e ~ 96(15 7 r) 1 /2L min 3 / 2 v ; 

If we scale this constraint using measured values of the constants, we obtain the relation 
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In our universe, black holes must have masses greater than about 6 x 10 13 g in order to 
last for iVbio = 10 33 atomic time scales, and must have masses less than about 2 x 10 14 
g in order to produce enough power (L min ). As a result, a biosphere could be powered 
by a black hole, although we have adopted somewhat optimistic requirements, e.g., the 
required luminosity is only L m i n , which is much less than a solar luminosity. The largest 
obstacle, however, is the production of black holes with this mass scale. 

Figure [6] shows the region of parameter space for which black holes can play the role 
of stars. To construct this diagram, we assume that black holes must live iVbio = 10 33 
atomic time scales and produce enough luminosity. For this latter requirement, we use 
the power intercepted from the Sun by the Earth (as a minimum value; dotted curve), 
the luminosity of a low-mass star (L ~ 1O~ 3 L ; dashed curve), and 1.0 L Q (solid curve), 
all scaled according to equation ( l57l) . If the black hole is required to have luminosity 
in the stellar range, then the allowed region of parameter space is highly constrained, 
in that the parameters (a, G) must have values quite far from those in our universe. In 
particular, the gravitational constant must be small (so that the luminosity is large), 
and the fine structure constant must also be small (so that atomic energy levels are 
low). If the necessary luminosity is determined by -L min0 = 10 17 erg/s, however, black 
holes can play the role of stars over a much wider range of parameter space. 



4-2. Degenerate Dark Matter Stars 

In principle, alternate universes can produce degenerate stars made of dark matter 
particles. Such stars could exist in our universe as well, although their formation is 
expected to be so highly suppressed that they play no significant role. This section 
considers the structure of these hypothetical objects in possible other universes. 

A degenerate star has the structure of an n = 3/2 polytrope, with the constant K 
in the equation of state given by 

K -^'^- < 60 > 

where is the mass of the dark matter particle. Since the constant K is specified, we 
can solve directly for the stellar properties. The mass-radius relation is given by 

M*Rl = ^o^W 8 G- 3 , (61) 

and the central density is given by 

32 G 3 m d s M* 
97rVo 2 h~ e 

For completeness, we note that the Chandrasekhar mass for these dark matter stars is 
given approximately by the expression 

f37rW 2 / ft c \ 3 / 2 

** = /*W-(gS->] »*. «») 

where fiQ ~ 2.714 for an n = 3/2 polytrope; this expression does not include 
general relativistic corrections (e.g., [T8]). For reference, note that a typical expected 



pc= :« • (62) 
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Figure 6. Allowed region of parameter space for the existence of black holes that can 
play the role of stars. The parameter space is the plane of the gravitational constant 
log 10 [G/Go] versus the fine structure constant log 10 [a/ao], where both quantities are 
scaled relative to the values in our universe. The allowed region lies under the curves, 
which are plotted here for three cases: The black hole luminosity is required to be 
greater than that of the Sun (solid curve), a low-mass star (dashed curve), and the 
solar luminosity intercepted by the Earth (dotted curve). The open triangle marks the 
location of our universe in this parameter space. 
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value for the dark matter particle mass, = lOOmp, implies that this mass scale 
M ch w 0.0007M o . 

For these stars, the luminosity is provided by annihilation of the dark matter 
particles. The annihilation rate per particle Ti is given by 

T 1 = n(cr d v) ps o- d ^n 4/3 /m d , (64) 

where <Jd is the cross section. To find the stellar luminosity due to dark matter 
annihilation, we must integrate over the star to find the total annihilation rate IV: 

F T = ^^/ 3 7o , where 7o = f* ?f/*d£ . (65) 
/io^d Jo 

For this n = 3/2 polytrope, 70 ~ 0.913. As a result, the total annihilation rate is 

given by Tp ps AjTd/3, where Nj- is the total number of particles in the star and Ti is 

evaluated at the stellar center. The corresponding stellar luminosity is then given by 

L -(I?) 4/3 ^^ G4M ' 1/V/3 - (66 » 

If the mass of the degenerate star were close to the Chandrasekhar mass, the luminosity 
would be enormous, and its lifetime would be short (see below). To put this in 
perspective, if we use reasonable values of the dark matter properties for our universe 
(ma = lOOmp and o"d = 10 -38 cm 2 ), then the mass required to produce L* = 1.0L© is 
only about M, ~ l(r 13 M Q ~ 10 20 g (about the mass of a large asteroid). As a result, for 
the range of parameter space for which these objects play the role of stars, the masses 
are far below the Chandrasekhar mass. 

If the dark matter star starts its evolution with initial mass M and later has a mass 
M(t) -C Mo, then its age At(M) is related to its current mass through the expression 

At(M) = § (^r ^5v*M- = i^t , (67) 

where is the luminosity of the star when it has mass M. For example, if M = 10 20 
g (the mass scale that generates = l.OL ), the time scale from equation fIBTj) is only 
about 100 days. In order for the time scale to be 1 Gyr, say, the mass scale must be 
about 3 x 10 16 g, and the corresponding luminosity is only ~ 10~ 13 L Q = 4 x 10 20 erg/s, 
i.e., still substantially larger than the expected value of L m ino- 

When the masses are well below the Chandrasekhar mass (see above), the star must 
satisfy two constraints. The first requirement is that the star is sufficiently luminous, 
which implies that 

„2_ ^4™ 25/3 

U = B ° ffdC ; 7 md M 11/3 > L min0 (a/a ) 4 , (68) 

h 

where we have defined a dimensionless constant B, 
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Next we require that the stellar lifetime is sufficiently long. In rough terms, this 
constraint can be written in the form 

At(M) = — — M" 8/3 > hNbio (70) 

where we have not made the distinction between M and M* in using equation (1671) . The 
first constraint puts a lower limit on the mass, and the second constraint puts an upper 
limit on the mass. By requiring that both constraints be met simultaneously, the mass 
can be eliminated and a global constraint can be derived: 

. ^ (<wy% Ca pq^y* (^y* (71) 

v a / \ he J L min0 \a^mdh / \ N bioJ 

where the constant Cb = (3/8) 11//8 .B~ 3 / 8 ~ 1.75. This result defines the parameters 
necessary for dark matter stars to play the role of ordinary stars (keep in mind that the 
formation of these bodies remains a formidable obstacle). The luminosity is determined 
by the dark matter annihilation cross section, which is independent of the constants 
that determine the physical structure of the star. As a result, the parameter space of 
constants (a, G) considered here always contains a region where these stars can operate: 
For fixed properties of the dark matter (md and ad), equation ( 1711) delineates the portion 
of the a — G plane that allows these degenerate dark matter objects to act as stars. 
On the other hand, one can use equation (I7T1) to constrain the allowed dark matter 
properties for given values of a and G. 



4-3. Other Possibilities for Unconventional Stars 

If the nuclear burning constant C = 0, then baryonic objects can still, in principle, 
generate energy in a variety of ways. In the absence of nuclear reactions, stellar 
bodies will often tend to form degenerate configurations, analogous to white dwarfs 
in our universe (provided that their mass is below the relevant Chandrasehkar mass 
scale). These degenerate objects can generate energy through several channels, including 
residual heat left over from formation, proton decay, and dark matter capture and 
annihilation. 

In the latter case, dark matter particles are captured by scattering off nuclei (which 
could be simply protons in a universe with no nuclear reactions). After a scattering 
event, the recoil energy of the dark matter particle can be less than the escape speed 
of the star, and the particle can be captured. After capture, the dark matter particles 
sink to the stellar center, where they collect until their population is dense enough for 
annihilation to balance the incoming supply of particles. The star thus reaches a steady 
state, where the luminosity is given by the total capture rate. This process has been 
discussed previously in a variety of context, including as a solution to the solar neutrino 
problem [19] and as a means to keep white dwarfs hot beyond their cooling times [20J. 

The capture rate of dark matter particles is given by 

r = Udma^drnVrel , (72) 
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where n^m is the number density of dark matter particles, cr^m is the total cross section 
for capture subtended by the star, and v re [ is the relative velocity These quantities 
depend on dynamical structure (distributions of density, velocity, angular momentum) 
of the background halo of dark matter |21j . In our universe, for example, the capture 
rate of dark matter particles by white dwarfs is of order T ~ 10 25 s~ x [20J. With the 
capture rate specified, the corresponding luminosity is given by 

L* = f v m d Y , (73) 

where is the mass of the dark matter particles and where the efficiency factor /„ 
takes into account energy loss from the star due to some fraction of the annihilation 
products being neutrinos. 

In this scenario, the luminosity depends on the number density of dark matter 
particles in the background (in the galactic halo in the context of white dwarfs in our 
universe). This density is independent of stellar properties. In a similar vein, the time 
scale over which the luminosity can be maintained depends on the overall supply of dark 
matter particles; this quantity is also independent of stellar properties. Thus, for any 
values of the constants (a,G), considered here as the relevant parameters that specify 
stellar properties, a universe can have the proper values of dark matter densities and 
cross sections so that degenerate stars can serve in place of nuclear burning stars. The 
specification of the allowed parameter space depends on more global properties of the 
universe, however, and is beyond the scope of this paper. 

5. Conclusion 

In this paper, we have developed a simple stellar structure model (§2) to explore the 
possibility that stars can exist in universes with different values for the fundamental 
parameters that determine stellar properties. This paper focuses on the parameter 
space given by the variables (G, a, C), i.e., the gravitational constant, the fine structure 
constant, and a composite parameter that determines nuclear fusion rates. The main 
result of this work is a determination of the region of this parameter space for which 
bona fide stars can exist (§3). Roughly one fourth of this parameter space allows 
for the existence of "ordinary" stars (see Figure [5]). In this sense, we conclude that 
universes with stars are not especially rare (contrary to previous claims), even if the 
fundamental constants can vary substantially in other regions of space-time (e.g., other 
pocket universes in the multiverse). Another way to view this result is to note that the 
variables (G, a, C) can vary by orders of magnitude from their measured values and still 
allow for the existence of stars. 

For universes where no nuclear reactions are possible, we have shown that 
unconventional stellar objects can fill the role played by stars in our universe, i.e., 
the role of generating energy (§4). For example, if the gravitational constant G and the 
fine structure constant a are smaller than their usual values, black holes can provide 
viable energy sources (Figure [6]). In fact, all universes can support the existence of stars, 
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provided that the definition of a star is interpreted broadly. For example, degenerate 
stellar objects, such as white dwarfs and neutron stars, are supported by degeneracy 
pressure, which requires only that quantum mechanics is operational. Although such 
stars do not experience thermonuclear fusion, they often have energy sources, including 
dark matter capture and annihilation, residual cooling, pycnonuclear reactions, and 
proton decay. Dark matter particles can also (in principle) form degenerate stellar 
objects (see §4). 

In order to assess the suitability of non-nuclear power sources, one must specify how 
much power is required, and for how long. In this work we have used the power that 
Earth intercepts from the Sun as the minimum benchmark value L m i n0 , and scaled the 
necessary power according to equation (I5T1) to account for variations in the fine structure 
constant; similarly, the required amount of time is taken to to ~ 1 Gyr, scaled by the 
atomic time of equation ( T4"9l) . These choices are not definitive, and hence alternative 
scalings can be explored. 

The issue of alternate values for the fundamental constants, as considered herein, 
is related to the issue of time variations in the constants in our universe. However, 
current experiments place rather strong limits on smooth time variations, with time 
scales exceeding the current age of the universe (see the review of |22j). Another 
possibility is for the constants to have different values at other spatial locations within 
our universe, although this scenario is also highly constrained |23j . 

This paper has focused on stellar structure properties. An important related 
question (beyond the scope of this work) is whether or not stellar bodies can be readily 
made in universes with varying values of the constants. Even if the laws of physics allow 
for stellar objects to exist and actively burn nuclear fuel, there is no guarantee that such 
bodies will be produced in significant numbers. In our universe, for example, there is 
a moderate mismatch between the mass range of possible stars and the distribution 
of masses of stellar bodies produced by the star formation process. At the present 
cosmological epoch, star formation produces objects over the entire possible range of 
stellar masses, with additional bodies produced in the substellar range (brown dwarfs). 
The matching is relatively good, in that the fraction of bodies in the brown dwarf range 
is small, only about 1 out of 5 [24] . Since the masses of these objects are small, the 
fraction of the total mass locked up in the smallest bodies is even smaller, less than 5 
percent. On the other hand, nearly all of the stars in our universe have small masses. 
As one benchmark, only about 3 or 4 out of a thousand stars are larger than the 
~ 8M threshold required for stars to experience a supernova explosion, whereas stellar 
masses can extend up to ~ 1OOM . The high mass end of the possible mass range is 
thus sparsely populated. The corresponding match between the range of allowed stellar 
masses and the mass range of objects produced can be quite different in other universes. 

In future work, another issue to be considered is coupling the effects of alternate 
values of the fundamental constants to the cosmic expansion, big bang nucleosynthesis, 
and structure formation. Each of these issues should be explored in the same level of 
detail as stellar structure is studied in this work. With the resulting understanding of 
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these processes, the coupling between them should then be determined. 

Finally, we note that this paper has focused on the question of whether or not stars 
can exist in universe with alternate values of the relevant parameters. An important and 
more global question is whether or not these universes could also support life of some 
kind. Of course, such questions are made difficult by our current lack of an a priori 
theory of life. Nonetheless, some basic requirements can be identified (with reasonable 
certainty): In addition to energy sources (provided by stars), there will be additional 
constraints to provide the right mix of chemical elements (e.g., carbon in our universe) 
and a universal solvent (e.g., water). These additional requirements will place additional 
constraints on the allowed region(s) of parameter space. 
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